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Abstract 

The perturbation method in supersymmetric quantum mechanics (SUSYQM) is used 
to study the spheroidal wave functions' recurrence relations, which are revealed by the 
shape-invariance property of the super-potential. The super-potential is expanded by the 
parameter a and could be gotten by approximation method. Up to the first order, it has 
the shape-invariance property and the excited spheroidal wave functions are gotten. Also, 
all the first term eigenfunctions obtained are in closed form. They are advantageous to 
investigating for involved physical problems of spheroidal wave function. 
PACs:11.30Pb; 04.25Nx; 04.70-s 

1 Introduction of the spheroidal functions 

The spheroidal wave equations are extension of the ordinary spherical wave equations. There 
are many fields where spheroidal functions play important roles just as the spherical functions 
do. So far, in comparison to simpler spherical special functions (the associated Lengdre's 
functions ) their properties still are difficult for study than their counterpart 
Their differential equations are 



d_ 
dx 



2^ d 



dx 



+ E + ax^ 



m 



1 



x^ 



e = Ox G (-i,+i). 



They only have one more term aX^ than the spherical ones (the associated Lendgre's equa- 
tions). 

This is a kind of the Sturm-Liouville eigenvalue problem with the natural conditions 
that Q is finite at the boundaries x = ±1. The parameter E can only takes the values 
£^0) El, ... , En, . . ., which are called the eigenvalues of the problem, and the corresponding 
solutions (the eigenfunctions) Oq, 6i, . . . , 6nj ■ ■ ■ are called the spheroidal wave functions 

m-m- 
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Under the condition a = 0, they reduce to the Spherical equation and the solutions 
to the Sturm-Liouville eigenvalue problem are the associated Legendre- functions P^(x)(the 
spherical functions) with the eigenvalues En = n(n + l), n = m + 1, m + 2, . . .. Though the 
spheroidal wave equations are extension of the ordinary spherical wave functions equations, 
the difference between this two kinds of wave funtions are far greater than their similar ity[T]. 

1. The spherical wave functions have many good property such as 

(a) the Legendre-functions Pn{x) = Pni^) are polynomials; 

(b) Back to the variable 9 with x = cos 9, the recursion relation among the Legendre- 



functions could be written as: 
Pn(cos( 



-ncos( 



smf 



d 
d9 



P„_i(cos( 



n 



1,2,3,... , 



(2) 



so, all Pn{cos9) could be deduced from the first or ground function Pq from the 
recursion relation 

(c) all the associated Legendre-functions P^{cos9) could be derived from by 

„ . „ d 



Pn (cos 9) 



-n cos t 



sml 



d9 



Pn-i{cos9), n = m + l,m + 2,m + 3, (3) 



2. On the contrary, very little are known about the properties of the spheroidal wave 
functions: 

(a) All spheroidal wave functions can not be polynomials at all. 

(b) Whether or not there are the similar recurrence relations between them is not 
clear, at least they have not been found upon to now. 



The recursion relations of the spherical wave functions, the associated Legendre func- 
tions reflect the the factorization properties of the Legendre equations [I6]. This shows that 
the Legendre equations are solvable. The Super-symmetry quantum mechanics tells that the 
roots of the solvable properties of differential equations are the super-invariance of their super- 
potential. In this paper, we testily investigate whether the spheroidal differential equations 
have the similar shape-invariance property. On accounting of the fact that their differential 
equations are difficult to treat, we only can rely on the approximate method of small param- 
eter a. The results may make one happy: they have shape-invariance property upon to the 
first term approximation. 



2 Brief review of the theory of SUSY in solving eigenvalue- 
problem 

In recent years, supersymmetric quantum mechanics have attracted tremendous attention 
for solvable potential problems. They not only provide clear insight into the factorization 
method of Infeld and Hull [16] , but also make great improvement in solving the differential 
equations. See reference [2] for review on its development. 
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In supersymmetric quantum mechanics [Il]-[l5], the Schrodinger equation is 

H-^l^- = -^y^ + V-{x)i;-{x) = E^^- (4) 

with % = 2m = 1. Here one considers the case of unbroken supersymmetry which demands 
that the ground state is nodeless with zero energy, and supposes that the superpotential 
W{x) is continuous and differentiable. The super-potential W{x) satisfies the equation 

V-{x)=W'^{x)-W'{x). (5) 

With the introducing the two operators 

A=^ + W{x) , = -^ + W{x), (6) 

the corresponding Hamiltonian H~ have a factorized form 

H- = A^A . (7) 
The partner potential V^{x) of is determined by the superpotential W{x) as 

V+{x) = W^{x) + W'{x), (8) 
and the corresponding Hamiltonian 

H+ij;+ = — + V- {x)7p+ (x) = E+iP+ (9) 

has a factorized form 

H+ = AA^ . (10) 

The Hamiltonians H'^ and have exactly the same eigenvalues (or the same spectrum) 
except that has an additional zero energy and the related eigenstate, that is, 

£;J-)=0, e\^\=E^-\ ^i+\cx^^(-) , ^t^{+)ocV^(-\ , n = l,2,.... (11) 

The pair of SUSY partner potentials V^{x) are called shape invariant if they are similar in 
shape and differ only in the parameters, that is 

V+{x;ai) = V-{x;a2) + R{ai), (12) 

where ai is a set of parameters, a2 is a function of ai (say a2 = /(ai)) and the remainder 
R{ai) is independent of x. One can use the property of shape invariance to obtain the ana- 
lytic determination of energy eigenvalues and eigenfunctions [H]-(15j. Thus for an unbroken 
supersymmetry, the eigenstates of the potential V~ (x) are: 

n 

E^=0, E- = Y,R{ak) (13) 

k=l 



oc exp 



(14) 



W{y,ai)dy 

^-=A^{x,ai)^>-_^{x,a2), n= 1,2,3,... (15) 
Therefore, the shape invariance condition actually is an integrability condition. 
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3 the eigenvalues and eigenfunctions of the excited states in 
first order 



In the following, we rewrite the equations ([T]) in the Schrodinger form. Though the form ([T]) 
is more familiar for research, the problem is easier to solve in the original differential equation 
than in the equation ([1]) . The original form is obtained from the eq. ([1]) by the transformation 



cos 9, 



that is. 



1 d 



+ a cos 



m 



sm 



e 



-EG 



(16) 



(17) 



the corresponding boundary conditions become is finite at 9 = 0, vr. 
Prom eq. ()17p and by the transformation 



G 



we could get the Schrodinger form as 



d^^ 



+ 



- + a cos 
4 



sin 2 9 



m 



sm 



and the boundary conditions 



The equations (|19p show the potential is 

1 



0. 



^' = 



V{9, a, m) 

The super-potential W could be determined by 



- — a cos 
4 



+ 



m 



sm 



W' = V{9, a, m) - Eq 



1 m 

■T + — 



4- — a cos^ 9 — En 



sm 



(18) 



(19) 



(20) 



(21) 



(22) 



Eq is the the first eigenvalue (the ground state energy) of the equations ()19p . Subtracting 
£"0 to make it possible to factorize eqns. (|19p . Actually, it is difficult to find the solutions of 
the equation ()22p . So the perturbation methods come naturally, which have been detailed in 
the reference [18j where the perturbation method in supersymmetry quantum has been used 
to resolve the ground state function of the spheroidal equations in approximation of little 
parameter a. Here is brief accounting of the results. One just expands the super-potential 
W , the central concept of supersymmetry quantum, in the series form of the parameter a. 
Because the potential is already in the series of a, one could solve the series forms of the 
super-potential W by solving term by term. Also note the fact the ground energy Eq must 
be in the series form of the parameter a and could be obtained by the boundary conditions. 
So, the super-potential W and Eq could be expanded as series of the parameter a, that is, 



W = Wq + aWi + a^W2 + a^W^ + 



(23) 
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n=0 



(24) 



The perturbation equation becomes 



n=0 



1 m 

■T + — 



sm 



oo 

acos"^ e -J^^ona'' (25) 

n=0 



There are two lower indices in the parameter E^n with the index refereing to the ground 
state and the other index n meaning the nth term in parameter a. The last term J2'n=o '^^^onOt^ 
is subtracted from the above equation in order to make the ground state energy actually zero 
for the application of the theory of SUSYQM. Later, one must add the term to our calculated 
eigen-energy. 

In the reference [T8], Wq, Wi, W2 are 



Wo 
Wi 

W2 



m + ^ ) cot 0, £^00 



sm y cos u 
2m + 3 

— sin 9 cos / 



+ 



m(m + 1). 



sm a cos ( 



(2m + 3)3(2m + 5) (2m + 3)2(2m + 5) 
The ground eigenfunction upon to the second order becomes 



exp 
A^'exp 

{sinOy 



WdO 



J WodO - a J Wide - j 



dO 



exp 0{a^ 



(26) 
(27) 
(28) 

(29) 
(30) 



exp 



a sm y 
4m + 6 



exp 



sin^ ( 



sin^ I 



2(2m + 3)3(2m + 5) 4(2m + 3)2(2m + 5) 



exp 0{a^ 



(31) 



Because the good result of the quantity Wi in the form of the equation (p7|) , the potential 
upon to the first order have the property of the shape-invariance. In order to see this property 
clearly, we rewrite the super-potential W as 



W = AiWq + aBiWi + 0{a^ 



(32) 



then 



V-{AuBi, 



W^ - W' 
{m + \fAl-{m + ]^)Ai 



+ 



V+{AuBi, 



Bia 
2m + 3 " 
W^ + W' 



(m + lfAl 



+ 0{a^) 



Bia 
2m + 3 



[(2m + l)Ai + 2] cos^ 6 



(33) 
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+ 



2m + 3 ^ 2-^ \ 
V.{A2,B2,e) + R{Ai,Bi)) 

{m + ^fAl-{m + ^)A2 



+ 0(a2) 



Bia 
2m + 3 



B2a 
2m + 3 



[(2m + l)Ai - 2] cos^ 9 



(34) 



[(2m + 1)^2 + 2] cos^ 6* 



+ 



B2a 
2m + 3 



1 



m 



+ 0(0^ 



+i?(^i,5i) +0(a^ 



where 



(m + 



^2 + 5 



2m + 3 



-a 



with ^1 = i?i = 1. From the above equations, it is easy to get 



(m+^)^2 

B2 

R{AuBi) 



/i(^i) = (m + i)Ai + l 



h{Ai,A2 



(2m + 1)^1 - 2 
(2m + 1)^1 +4 



Si 



Define 



fc+i 



(m + -Mi + l--(i?2+i?i) 



/i(^fc) =^fc + 



m + i 



R{Ak, Bk, Bk+i) 



(2m + + 1 - f (^fc+i + ^fc) 



then the n+lth energy (or the nth excited energy)upon to the first order is 

n 

En = ^RiAk,Bk,Bk+i) 



k=l 

n 



J2 ((2m + l)Ak 



a 



k=l 



n(2m + n + 1) 



2m + 3 
a 



{Bk+i + Bk) 



2m + 3 



[2 (5i + ^2 + . . . + Bn) -Bi + B, 



n+lj 



By the relation of eq. (|^T|) , we get 

{Bk+i + Bk) 

therefore, 



2k 



2m + 3 



{Bk — Bk+i) 



{Bi + B2 + ... + B„ 



1 



[(2m + 3)^1 - (2n + 2m + 3)5„+i] . 



(35) 

(36) 

(37) 
(38) 
(39) 

(40) 

(41) 
(42) 

(43) 

(44) 
(45) 

(46) 
(47) 
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We could get Bn+i from the relation ([^T|) . that is, 

(2m - l)(2m + l)(2m + 3) 



B. 



so, 



n+l 



(2n + 2m - 1) (2n + 2m + 1) (2n + 2m + 3) 



5i 



(48) 



2(2m + 3) 
n(2m + n + 1 
a 



[(2m + l)5i - (2n + 2m + 1) + 0(0^ 



2(2m + 3) 



_ (2m-l)(2„.+ l)(2„. + 3) 

^ ' (2n + 2m - 1) (2n + 2m + 3) 



0{a^). (49) 



As stated in the reference |18] . before, the last term in eq. ()25p must be added to -E„ for 
actual calculation, that is. 



+ ii^oo + Eoia = n(2m + n + 1) + m(m + 1) 



a 



2m + 3 



a 



2(2m + 3) 



(2,„-l)(2,„ + l)(2,„ + 3) 

^ ^ (2n + 2m - 1) (2n + 2m + 3) 



a 



n{2m + n + 1) + m(m + 1) — — 



1 



(2m - l)(2m + 1) 



(2n + 2m - 1) (2n + 2m + 3) 



(50) 
(51) 



Notice that our formula is different from that of the reference [T] , the relation is that our 
m + n is equivalent to / in the reference [Tj, that is, n = Z — m and an overall negative sign 
between them. Then we can compare the results, which is same: 



+ £^00 + Eqiu + 0(a) = {I - m){l + m + 1) + m(m + 1) 



a 



1 



(2m - l)(2m + 1) 
{21 - 1) {21 + 3) 



/(/ + !) 



a 



(2m - l)(2m + 1) 



{21 - 1) (2/ + 3) 



The nth excited state upon to the first order is obtained by 
i^-{9,A^,B^)=AH-_^{e,A2,B2) 



oc 



A\e, Ai, Bi)A\e, A2, B2)... A\9, An, B„)V'o"(^, ^„+i, Sn+i) 



(52) 



(53) 
(54) 



with 



%{e,An+uBn+i) = (^m^)(™+^)^n+i 



exp 



aBn+i sin^ 
4m + 6 



*expO(a2), (55) 



and 



a 



A^{e, An,Bn) = -4 - ("1 + \)K cot e - —^n 

ad 2 4m + b 



Br, sin 26 



(56) 



where An, An+i satisfy the relation of ([in|) . ([^T]) . 
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4 conclusion and discussion 



In the paper, the first order term of the super-potential obtained has the good property of 
shape-invariance by the formula of ()27p . this in turn makes it easy to get the pleasant results 
of the excited energy and the corresponding excited state functions from the counterparts of 
the ground state. The results are much more like case of the associated-Lengdre functions. 
They are new and interesting. The method can also be used to solve higher order terms of 
the super-potential and the related super-invariance problem of the spin-weighted spheroidal 
functions. The methods apply to the spin- weighted spheroidal functions for the case s ^ 0. 
These cases have been studied too and will be reported. 
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